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Abstract 

We develop an approach for the treatment of one-dimensional 
bounded quantum-mechanical models by straightforward modifica- 
tion of a successful method for unbounded ones. We apply the new 
approach to a simple example and show that it provides solutions to 
both the bounded and unbounded type of models simultaneously. 

1 Introduction 

The Riccati-Pade method (RPM) yields accurate eigenvalues and eigenfunc- 
tions of separable quantum-mechanical models [1-13]. The approach is based 
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on a rational approximation to a modified logarithmic derivative of the eigen- 
function and the best fit occurs when the eigenvalue is a root of a Hankel 
determinant [1-13]. 

The roots of the Hankel determinant approach the energies of bound 
states and resonances as the determinant dimension increases. The RPM 
does not require explicit specification of the boundary condition, the ap- 
proach commonly selects the physical one automatically for each problem. 
The resulting eigenvalues always correspond to the correct asymptotic be- 
haviour at infinity [1-11]. As far as we know the RPM is the only approach 
that applies exactly in the same form to both bound states and resonances. 

If the potential-energy function has poles at two points, then the roots 
of the Hankel determinant approach the eigenvalues of the problem with 
Dirichlet boundary conditions at such coordinate values [8,12,13]. We call 
these boundary conditions "natural" . 

In some cases one is interested in that the eigenvalue equation satisfies 
"artificial" bounday conditions. For that reason, in this paper we propose a 
modification of the RPM to treat Dirichlet boundary conditions at arbitrary 
coordinate locations. In Sec. [2] we introduce the RPM for one-dimensional 
models and suggest how to force the desired boundary conditions. In par- 
ticular we concentrate on a linear potential that has proved useful for the 
treatment of some physical problems. In Sec. [3] we show results for the cho- 
sen eigenvalue equation, and in Sec. H] we summarize the main features of 
the RPM and draw some conclusions. 
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2 The method 



We introduce the RPM by means of a simple one-dimensional problem of 
the form 

Y"(x) + Q(x)Y(x) = 0, Y(0) = (1) 

that depends on an adjustable parameter that is necessary to satisfy the 
other boundary condition which we will specify below. For example, in the 
case of a dimensionless Schrodinger-like equation Q(x) = V(x) — E, where 
the energy E is the adjustable parameter. For concreteness we restrict to 
this case from now on. 

In order to apply the RPM we define the modified logarithmic derivative 

where the function g(x) is chosen so that f(x) is analytic at x = and 
therefore can be expanded in a Taylor series 

oo 

/(*) = E If' (3) 

j=0 

Notice that the coefficients fj depend on E. The RPM is based on the trans- 
formation of the power series ([3]) into a rational function or Pade approximant 
that satisfies 

V M n r3 M+7V+1 

[M/N}(x) = |£Li- = £ ./>H • ^ (4) 

2^=0 °j x j=o 

where M = N + d, d = 0, 1, . . .. Notice that the rational ansatz has just 
M + N + 1 adjustable parameters dj and bj to fit the first M + N + 2 
coefficients of the Taylor series ([3]). The additional requirement determines 



3 



the value of E clS cL root of the Hankel determinant: 

HUE) = \fr +3+d+1 (E)\ hj=OAN = 0, (5) 

where D = N + 1 is the dimension of the Hankel matrix. Each Hankel 
determinant is a polynomial function of E and we expect that there is a 
sequence of roots E^ D,d \ D = 2, 3, . . . that converges towards the value of E 
consistent with the second boundary condition. 

Commonly, the Hankel quantization condition (jSj) provides the eigenval- 
ues consistent with the bound states (Y(x — > oo) = 0) or the resonances 
embeded in the continuum (outgoing or incoming waves). The RPM au- 
tomatically selects the eigenvalues that are consistent with such "natural" 
boundary conditions [1-11]. 

If the potential-energy function exhibits poles, then the RPM automat- 
ically selects Dirichlet boundary conditions at the corresponding coordinate 
points. For example, when V(x) = Vosec(x) 2 the RPM selects the boundary 
conditions F(±tt/2) = [8], and Y(±R) = when V{x) = a 2 x 2 /(l-x 2 / R 2 ) 2 
[12,13]. 

In some cases one wants to force boundary conditions that are not re- 
lated to singular points in the potential-energy function. Suppose that we 
are interested in the differential equation ([T]) with the boundary conditions 
Y(0) = Y(l) = 0. We can force such "artificial" boundary conditions by 
means of a properly chosen function g(x) in equation ([2]). In fact, the func- 
tion g(x) = x(l — x) introduces poles at x = and x = 1 into the differential 
equation for f(x) that we can rewrite as 

x(l - x)f'(x) + 2(1 - 2x)f(x) - x(l - x)f(x) 2 - x{l - x)Q(x) + 2 = (6) 
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In this way we expect to obtain the eigenvalues consistent with those bound- 
ary conditions. 

For simplicity we consider 

Q(x) =e-Xx (7) 

A motivation for this choice is that the resulting differential equation and 
boundary conditions are related to a simple model for the study of electrons 
in a crystal under the effect of an electric field [14]. The Schrodinger equation 

- + eFX$(X) = E${X), $(0) = = (8) 

provides the states and energy levels of an electron of mass m and charge e 
in a box of impenetrable walls at X = and X = L (that mimics the finite 
size of the crystal) under de effect of an electric field of strength F [14]. This 
extremely simple model has also been useful in the study of the tail of the 
density of states of a disordered system in the presence of an electric field [15]. 
By means of the change of variables X = Lx and <&(Lx) = Y(x) one obtains 
the differential equation flTJ with the coefficient ([7]) where A = 2mL 3 Fe/h 2 , 
and e = 2mL 2 E/h 2 . 

Another reason for the choice of such example is that one can write its 
solutions exactly in terms of the Airy functions Ai(z) and Bi(z): 

Y(x) = N [Bi{-e)Ai{\x - e) - Ai(-e)Bi(Xx - ej] (9) 

where iV is a normalization factor, and the dimensionless eigenvalues e n , 
n — 0, 1, . . . are given by the quantization condition 

Bi(-e)Ai(X - e) - Ai(-e)Bi(X - e) = (10) 
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3 Results 



The application of the RPM is straightforward: we obtain as many coeffi- 
cients fj(e) as necessary from the differential equation for f(x), construct the 
Hankel determinants H^(e) D = 2, 3, . . . and calculate their roots. We expect 
these roots to converge towards the eigenvalues of the differential equation 
with the boundary conditions mentioned above. 

Table [JJ shows sequences of roots of the Hankel determinants that already 
converge towards the exact eigenvalues given by equation ({TO]) when A = 1. 
As in previous applications of the RPM we appreciate that the rate of con- 
vergence of the Hankel sequences decreases as the energy increases because 
the denominator of the rational approximation (J3J) requires greater values of 
iV to accomodate the increasing number of zeros of the solution Y(x). We 
clearly see that this modification of the RPM enables one to solve eigenvalue 
equations with "artificial" Dirichlet boundary conditions. 

The Hankel determinants are polynomial functions of the eigenvalues and 
display many more roots than those that we choose to build the sequences 
that converge towards the actual eigenvalues of the given problem. One of 
the features of the RPM is that an increasing number of roots cluster around 
the eigenvalues as D increases. For the simple example chosen here there are 
only two roots that approach a given eigenvalue as D increases ( at least for 
D < 16). Fig. [TJ shows log \6q PPTOX (D) — e^° act \ for these two sequences. 

In the present case the Hankel determinants exhibit other roots than those 
mentioned above. They correspond to the "natural" boundary condition 
Y(x — > oo) = with eigenvalues given exactly by the quantization condition 
Ai(— e) = 0. The choice of g(x) suggests that we are looking for a solution of 
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the form Y(x) = x(l — x)e~ / ^ dx but the RPM also selects a solution of the 
form Y(x) = xe~ / dx with the "natural" boundary condition at infinity. 
The rational approximation to f(x) — f(x) + 1/(1 — x) absorbs and removes 
the pole at x — 1 and produces sequences of roots that converge towards the 
solutions of the unbounded problem (0 < x < oo). Table [2] shows some of 
these eigenvalues for A = 1. Curiously, more roots cluster around a given 
eigenvalue of the unbounded model than of the bounded one. Fig. [2] shows 
log \eQ Pprox (D) — €Q Xact \ for all the sequences that appear when D < 16. 

The function g(x) = x is more convenient for the "natural" boundary 
conditions and, consequently, the sequences of roots of the Hankel determi- 
nants exhibit greater convergence rate. Fig. [3] shows log \€q PPTOX (D) — ^ xact \ 
for the optimal sequences for both choices of g(x). When g(x) = x(l — x) 
the rational approximation to f(x) has to remove the wrong zero at x = 1 
and, for this reason, the rate of convergence of the RPM is slightly smaller. 
Notice that the solution Y(x) that satisfies Y{x) = Y(x — > oo) = does not 
have a cero at x = 1. 

4 Discussion 

Simple models of bounded quantum-mechanical systems have proved useful 
for the study of several physical phenomena [16] (and references therein). 
The modification to the RPM proposed here is suitable for bounding a sys- 
tem between impenetrable walls that force Dirichlet boundary conditions at 
their locations. The numerical results of the preceding section show that 
the convergence rate of the modified RPM is as remarkable as in the case 
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of the unbounded and naturally bounded systems [1-13]. A curious fea- 
ture of present application of the RPM to a bounded model is that the 
approach also provides the eigenvalues of the unbounded one. This outcome 
is a consequence of the fact that the RPM automatically selects the correct 
asymptotic behaviour at infinity of the solution to the differential equation. 
In all the cases studied that asymptotic behaviour coincided with the one 
required by physical reasons (vanishing at infinity, incoming or outcoming 
waves, etc) [1-11]. 

The transformation of the Schrodinger equation into a Riccati one has 
proved suitable for the application of the quasilinearization method (QLM) 
to quantum mechanics [17-22]. Regarding the calculation of the resonances 
of a quartic anharmonic oscillator the RPM [6] proves to be more accurate 
than the QLM [22]. 

The main ideas behind the RPM have recently proved useful for the 
treatment of two-point nonlinear equations [23] of interest in some fields of 
physics [24-26] . The resulting approach called Hankel-Pade method (HPM) 
appears to be an alternative accurate tool for the determination of unknown 
parameters of the theory that are consitent with the desired asymptotic be- 
haviour of the solution of the nonlinear differential equation [23]. 
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Table 1: First four eigenvalues of the bounded model 



D 


eo 




2 


9 




3 


10.2 




4 


10.36 




5 


10.3679 


35 


6 


10.36848 


39.3 


7 


10.368506 


39.89 


8 


10.36850713 


39.97 


9 


10.368507161 


39.978 


1 n 


1 n ^£8^(171 fi1 897 
lU.OOoOLw lOloZi 


qq Q787 

oy.y ( 1 


11 


10 3685071618362 


39 97874 


12 


10 368507161836336 


39 9787445 

t7>t7 1 ^ 1 Jl "_7 


13 


10 3685071618363371 


39 97874477 


14 


10 368507161836337126 


39 9787447892 

U t7 • t7 1 W 1 JT | Ui7 4J 


15 


10 368507161836337127 


39 97874478986 


16 


10 368507161836337127 


39 978744789882 


Exact 


in QdQ rn7i £ 1 oq^qq ti 07 
IU.0000U / 10100000 f iz f 
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D 




^3 


8 


81 




9 


88 




10 


89.1 




11 


89.3 


144 


12 


89.321 


156 


13 


89.3259 


157.9 


14 


89.3266 11 


158.31 


15 


89.326628 


158.39 


16 


89.3266340 


158.411 



Exact 



89.326634542478746080 



158.41378981431004871 



Table 2: First four eigenvalues of the unbounded model 



D 


eo 




4 


2.29 




5 


2.337 




6 


2.33808 


4.0 


7 


2.3381070 


4.083 


8 


2.33810740 


4.0878 


9 


2.3381074103 


4.087945 


i n 

1U 


9 Q'iSI 177/1 1 D/1 ^P. 




11 


2 3381 0741 045970 


4 D87949441 


12 


2 338107410459766 


4 0879494440 


13 


2 33810741045976702 


4.087949444129 


14 


2 3381074104597670382 

• - 7 -.7 v7 -L V 7 1 -L -L V 7 J- .7 ' 7 1 V7 1 V 7 -.7 V - 


4 08794944413093 

X • V 7 v - 1 ' 7 -L ' 7 A X J- -L - 7 V 7 7 -.7 


15 


2.3381074104597670385 


4.087949444130970 


16 


2 3381074104597670385 


4 08794944413097060 


Exact 


z.ooolU ( 41U4oy /Uooo 


/I I7S707 O/l 7 7 1 Q(707nf;i CP. 
4.Uo ( y4y4441oUy ( UOIOO 


D 






8 


5.1 




9 


5.50 




10 


5.520 




11 


5.52054 


6.74 


12 


5.5205591 


6.785 


13 


5.52055981 


6.7866 


14 


5.520559827 


6.786705 


15 


5.52055982808 


6.7867080 


16 


5.5205598280950 12 


6.786708086 



Exact 



5.5205598280955510591 6.7867080900717589988 




Figure 1: Sequences of roots of the Hankel determinants for the lowest eigen- 
value of the bounded model 




Figure 2: Sequences of roots of the Hankel determinants for the lowest eigen- 
value of the unbounded model 
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Figure 3: Sequences for the lowest eigenvalue of the unbounded model for 
g(x) — x(l — x) (circles) and g(x) = x (squares) 
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